Peculiar Velocity Anomaly from Forces Beyond Gravity? 
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We address recently reported anomalously large bulk flows on scales of 100/i -1 Mpc and beyond. 
These coherent motions of galaxies challenge the standard ACDM concordance model as well as 
a large class of competitive models of dark energy and modified gravity. If confirmed, they may 
support alternative models that include extra forces enhancing the growth of perturbations on large 
scales. In such scenarios, current observational constraints on structure formation favor an onset of 
the extra forces in recent times, as predicted, e.g., in growing neutrino models. For the illustration 
of the main effects, we employ intuitive phenomenological parameterizations. 



Recent observations of large-scale galaxy motions con- 
stitute one of the main challenges for the cosmological 
standard model In a Gaussian window of diameter 
100/i -1 Mpc, Watkins et al. Q find a coherent bulk mo- 
tion of 407 ± 81 km/s in conflict with the expectation of 
ps 200 km/s at the 2a level. Other independent results con- 
firm the presence of unexpectedly large bulk motions 0, 0| • 
Kashlinsky et al. Q investigate scales of « 300/i _1 Mpc, 
where the expectation is even lower, obtaining the drastic 
result of 600-1000 km/s. Despite of the large uncertainties 
still present today, such values have the potential of form- 
ing a highly significant anomaly for the ACDM model in 
the future. 

Yet, the standard ACDM model has passed a series of 
stringent tests mainly from measurements of the cosmic mi- 
crowave background (CMB), from galaxy, Lyman-a, weak 
lensing, and supernova la surveys [5j, [y, |7j, [8J, l9|. It as- 
sumes a spatially fiat universe essentially made up from 
dark energy in the form of a cosmological constant A, 
cold dark matter (CDM), baryons and radiation. Accord- 
ing to the standard picture, matter perturbations from a 
nearly scale-invariant primordial spectrum grew solely due 
to Einstein gravity. Since the peculiar velocity field is in- 
timately connected to the growth of structure, a modi- 
fied growth history will typically affect the expected pe- 
culiar velocities. This can occur, e.g., in models of mod- 
ified gravity 13, brane- world models HI, 14. 15 1 
and models inclu ding extra couplings in the dark sector 

H E3, EE 0, 0, H2 0, 0] . 

This letter aims at illustrating that a modified growth 
history could achieve the large peculiar velocities suggested 
by the observations, while at the same time remaining 
consistent with further observational constraints. Some 
authors have already proposed explanations of the 
anomaly, namely supcrhorizon sized non-Gaussian, non- 
inflationary inhomogeneities 24|, cosmolo gica l-scale extra 
dimensions [l3| , and moving dark energy |25| . Instead of 
adopting any particular model, we parameterize modifica- 
tions to the standard model phenomenologically. Thereby, 
we point out the basic features an explanation of the 
peculiar velocity anomaly should have. We will also in few 
words discuss the idea of altering the primordial spectrum 
rather than the dynamics. It is our hope that this 
letter motivates future investigations based on concrete 
cosmological models. 



Before we introduce and study our phenomenological 
parameterizations in view of current observational con- 
straints, we collect the relevant ingredients from pertur- 
bation theory. Inhomogeneities in the metric induce devi- 
ations from the uniform Hubble flow. They are accounted 
for by the peculiar velocity field v(x). The bulk flow u is 
the average peculiar velocity in some volume defined by a 
window function W , 



u(x) 



dV(y)M/(x-y) 



(1) 



Every cosmological model predicts the mean square (u 2 ) 
for a window of given size and shape, which can be com- 
pared with observation. Throughout, we stick to statisti- 
cal homogeneity and isotropy, which for the Fourier trans- 
formed velocity field Vk implies 



(v k v k ,) = (2ir) 3 P v (k)S 3 (k-k'), 



(2) 



where we have introduced the peculiar velocity power spec- 
trum P v (k). It enables us to write 
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(3) 



for a Fourier transformed spherically symmetric window 
W{k). In linear perturbation theory, P v {k) can be obtained 
from numerical codes like CAMB [26| or CMBEASY [23 . We 
use the gauge-invariant velocity perturbation. 

It is instructive to relate the peculiar velocity power spec- 
trum P v (k) to the matter density power spectrum Ps(k). 
This is achieved with the aid of the continuity equation. In 
conformal Newtonian gauge, it reads 



-ik • v k + 30k, 



(4) 



where a dot denotes the derivative with respect to confor- 
mal time. We assume that the interactions of the matter 
fluctuations are dominated by gravity, with the metric per- 
turbation 0k as in [28j. In the Newtonian limit, 0k can 
be neglected. We define = dlogJk/dloga, such that 
Eq. (g]) reads 



-ik 



v k , 



(5) 



2 



with the physical Hubble parameter H. In linear perturba- 
tion theory, correlations between and 5\< are neglected. 
Then, performing the averaging of Eq. ([5]) and introducing 
the averaged growth factor /, 



/(*) 



1 dlogPf(fc) 

2 d log a 



yields the simple relation 



f 2 H 2 

P v (k)= J —P s {k). 



(6) 



(7) 



Equation tells us that larger bulk flows demand higher 
values of / or Pg(k). This, however, poses a serious obstacle 
for most cosmological models that reproduce the standard 
expansion history. On the one hand, once a model is cho- 
sen, the density power spectrum Pg{k) is constrained by 
various observations (like the CMB and galaxy surveys). 
Therefore, most cosmological models don't allow for dras- 
tic deviations from the ACDM power spectrum. On the 
other hand, for a large class of dark energy and modi- 
fied gravity models, / can be parameterized by f = fij^ 
with 7 constant in time . Linder and Cahn [30| showed 
that for models of uncoupled dark energy, 7 only slightly 
depends on the equation of state w, and that even when 
considering models of modified gravity, 7 typically varies 
at most w 20%, not enough to predict the observed bulk 
flows. Consequently, what at first was found as a challenge 
for the standard model ACDM, is in fact a problem for its 
most popular competitors as well. 

This should come as no surprise since the direct influ- 
ence of uncoupled dark energy, whether it be dynamical 
or a cosmological constant, is restricted to the evolution 
of the background. Similar expansion histories thus imply 
similar growth histories. This correspondence is absent in 
models with extra couplings enhancing the growth of mat- 
ter perturbations in recent times. For these models, the 
parameterization of / mentioned above in general becomes 
invalid. Going back to Eq. dH) and using Eqs. (0 and (|6|), 
we obtain: 



(« 2 ) = 



1 



87r 2 a 2 



(S) 



This formula shows what the continuity equation Q al- 
ready indicates: As P$ and P$ are evaluated at the present 
time, the bulk flows are sensitive to the growth of matter 
perturbations today. Hence, it is possible to reach large 
bulk flow expectations if there is significant current growth. 
This can be mediated by extra forces as present in models 
of coupled dark energy [B, E3, 0, O M, ED, [H, H3 • A 
rapid growth of perturbations would also contribute to the 
integrated Sachs- Wolfe effect (ISW), predicted too small 



for ACDM [3l|. If the late-time growth has set in very 



recently, the impact on the present power spectrum Pg 
may be moderate. There is thus a good chance to remain 
consistent with observational constraints on P$. 



In the following, we shall illustrate with the help of 
intuitive phenomenological parameterizations that ampli- 
fied late-time growth actually can succeed in resolving the 
peculiar velocity anomaly. We build upon the standard 
WMAP-5 best-fit ACDM model [§], keeping the back- 
ground evolution fixed, and introduce extra forces starting 
to act on the perturbations at some redshift 2*. We in- 
vestigate two different scenarios. First, an attractive extra 
force may directly act upon matter perturbations. Second, 
cold dark matter may only indirectly be affected by some 
additional gravitational potential caused by an extra force 
between other constituents (e.g. dark energy and neutri- 
nos, as in 17, 3, We will adjust the extra forces such 
that the observed bulk flows of [2j coincide with the theo- 
retical expectation and check whether we are still consis- 
tent with observational constraints for the cold dark matter 
power spectrum from the SDSS DR7 @|, the ISW [HI, (3!, 
and the CMB angular power spectrum from WMAP-5 33j . 
There are publicly available likelihood codes that can be 
implemented in COSMOMC [HI]. 

The first scenario is modeled by amplifying the coupling 
between matter perturbations. In the perturbation equa- 
tions, we modify the effective Newton constant G e ff: 



G 



off 



G for z > z* 
aG for z < z* 



(9) 



Here, G denotes the usual Newton constant and a > 1 the 
amplification factor. This can be interpreted as adding an 
extra force of strength a — 1 (compared to gravity) switched 
on at z*. Since the bulk flow anomaly is found on large 
scales, it is sufficient to restrict the modifications to scales 
k < fc*, where fc* is some large scale. Motivated by the bulk 
flow measurement by Watkins et al. 2] in a Gaussian win- 
dow of radius R = 50/i _1 Mpc, we chose fc* = tt/R. Zhao 
et al. 35] provide a CAMB variant for modified growth that 
makes the implementation of the modification ([9]) particu- 
larly simple. 

For varying values of z*, we chose a such that the bulk 
flow prediction from Eq. ([3]) equals the measured bulk flow 
u bs = 407±81 km/s. The results are shown in figure[T] A 
region of particular interest is around z* ~ 0.4 where dark 
energy starts to dominate in ACDM. There, at the latest, 
the altered dynamics of coupled dark energy models are 
expected to become important. Figure [T] suggests that the 
observed bulk flows could be due to an extra force, turned 
on at z* = 0.4 and two to four times stronger than gravity. 

We next confront the parameters favored by the bulk 
flow measurement with further observational constraints. 
The constraints from Lyman-a and weak lensing needn't 
be considered because we did not change small-scale dy- 
namics. We find that the WMAP-5 likelihood varies only 
slightly compared with the ACDM model for the parameter 
range considered here. But as can be seen in figure 2(a)| 
the SDSS DR7 imposes important constraints. All switch- 
on redshifts z* > 1 are excluded beyond 2a. These borders 
are sensitive to the choice of fc* that defines the scales at 
which the extra force is acting. Going to smaller fc* will 
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cate that both are actually larger than expected in the 
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FIG. 1: The parameter space for the parameterization ©. The 
shaded region marks the points consistent with the observation 
0| at < la. The continuous line highlights the points for which 
the measured value is hit exactly. 
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ACDM model. Ho et al. [31( use the relative amplitude A 



quantifying how strong the ISW was measured compared 
with the model-dependent expectation. Figure |2(b)| shows 
that the value of A improves in the parameterization §§§ 
and is closer to the ideal value A = 1 than it is for ACDM. 
The intuitive expectation is thus confirmed that current 
growth leads to a stronger ISW; the resulting values are 
even preferred by the observation. 

Let us now consider the second scenario, i.e., that matter 
is affected by an additional gravitational potential ex tra- 
Similar to the first scenario, we simply parameterize e xtra 
by 



for z > z* 
X for z < Z* 



(10) 



where x 1S a constant and z* & switch-on redshift. As be- 
fore, the modification is restricted to large scales k < fc*. 
The strength of the additional potential is most clearly 
quantified by the ratio x/4>i where (f> is the gravitational 
potential in ACDM evaluated at very large scales. Fig- 
ure[3]visualizes the parameter range for which the bulk flow 
measured by Watkins et al. 0] is reached. If the additional 



-e- 
x 




0.2 0.4 0.6 0.8 1 1.2 1.4 

Switch-on redshift z* 



FIG. 3: The parameter space for the parameterization IjlOj l. 
The region consistent with the observation of 0] is shaded, pa- 
rameters that predict the measured value exactly lie on the red 
line. 



FIG. 2: 



(a) The reduced \ 2 values from the SDSS DR7 matter 
power spectrum for the marked points in figure [T] The grey 
dashed lines mark the upper borders of the 1 and 2 a likelihood 
intervals respectively. The green dotted line refers to the ACDM 
best fit. (b) The measured ISW amplitude relative to the model 
prediction (uncertainties amount to « 25%). The green dotted 
line marks ACDM, the grey dashed line gives the ideal value 
A = l. 



considerably soften the constraints because current galaxy 
surveys lack sufficient statistics on very large scales (the 
largest scale entering the SDSS DR7 likelihood analysis [B| 
is k min = 0.02/i/Mpc). 

Besides the bulk flows, the ISW is a probe of the growth 
of matter perturbations. Remarkably, observations indi- 



potential is switched on at = 0.4 (equality of matter 
and dark energy), figure [3] indicates that xl4> should take 
values from approximately 1 to 3. 

The constraints from the SDSS DR7 power spectrum and 
the ISW are very similar to the ones already found for the 
first scenario and visualized in figures 2(a) and 2(b) Again, 
z* < 1 is strongly favored by the SDSS, and the ISW is in 
better agreement with observation than the ACDM pre- 
diction. The CMB, however, is notably affected. For the 
particular model with a sharp jump of O xtra, the predicted 
power spectra are excluded by the data at more than 2a. 
This issue depends, however, in a crucial way on the par- 
ticular time history and scale dependence of </> ex tra- One 
has to assess concrete models in order to see whether a 
scenario of this type can be consistent with WMAP-5 and 
nevertheless yield an enhanced bulk flow. 
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Finally, we briefly comment on the possibility of altering 
the primordial spectrum of perturbations in order to 
obtain the observed bulk flows without abandoning the 
standard dynamics. A significant shift of power on very 
large scales may in principle be motivated, e.g., by models 
of double inflation. Though an appealing idea at first 
sight, the CMB imposes stringent bounds. When the 
measured bulk flows are to be obtained without losing 
consistency with the CMB, extreme values have to be 
selected. Estimations suggest that one would have to go 
to scales far outside the horizon and to amplify the power 
spectrum by a large factor. Moreover, re gard ing the lack 
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of large-scale power in the CMB maps [3 
would seem more natural to assume less primordial power 
on the largest scales. 

In this letter, we have tried to point out which sorts of 
modifications to the standard cosmological model could be 
required in order to solve the peculiar velocity anomaly. 
Enhanced growth of matter perturbations in recent times 
z < 1 is a natural possibility. We have illustrated this with 
a rough parameterization of an extra force ([5]), which suc- 
ceeds in explaining the measured peculiar velocities while 
respecting further constraints. Bearing in mind the sim- 
plicity of our parameterization, it is to be expected that 
more elaborate cosmological models leading to enhanced 
growth can equally resolve the anomaly. We have shown 
that ISW observations even support the idea of enhanced 
growth; the predicted amplitude, too small in ACDM, 
comes closer to the observationally preferred value (see fig- 
ure 



2(b)| 



We have also considered the case in which the additional 
growth is not caused by an extra force directly acting on 
CDM particles, but rather indirectly via some extra grav- 
itational potential. This scenario as well can be in accor- 
dance with the peculiar velocity measurements and current 
constraints from large-scale structure. Our discussion has 
shown that scenarios of a fast-growing gravitational po- 
tential influence the CMB. It has to be studied carefully 
for every individual model whether the modifications are 
acceptable or not. 

If the significance of the peculiar velocity anomaly 
increases with future observations, this would be a strong 
sign for physics beyond the ACDM concordance model. In 
our examples, modifications to ACDM can only become 
effective in the recent epoch. This might suggest that 
the new physics is connected with the transition to dark 
energy domination. 

We thank Bjdrn Maltc Schafer for interesting and useful 
discussions and Gong-Bo Zhao for helpful comments on 
his CAMB variant. We acknowledge the use of the publicly 
available likelihood codes and data sets employed in this 
work. 
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